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We study the statistical properties of the eigenvalues of non-Hermitian operators assoicated with 
the dissipative complex systems. By considering the Gaussian ensembles of such operators, a hier- 
archical relation between the correlators is obtained. Further the eigenvalues are found to behave 
like particles moving on a complex plane under 2-body (inverse square) and 3-body interactions and 
there seems to underlie a deep connection and universality in the spectral behaviour of different 
complex systems. 



The random non-Hermitian operators play a significant role in the dynamics of variety of complex systems e.g. 
dissipative quantum systems Q, chaotic quantum scattering [^j, neural network dynamics ||, statistical mechanics 
of flux lines in superconductors with columnar disorder Q , classical diffusion in random media , biological growth 
problems || . The study of the statistical properties of their eigenvalues and eigenvectors therefore is relevant and has 
been of great interest in recent past. 

The random nature of a quantum operator reveals itself through a distribution of values for each of its matrix 
elements. The type of the distribution depends on the complexity of the system, resulting in a variety of non- 
Hermitian random matrices (NHRM) for various dissipative cases. In this letter, we develop a technique to study 
the spectral properties of the Gaussian ensembles of NHRM. The technique is based on a mapping of the eigenvalue 
distribution of a general NHRM to a non-stationary state of a classically integrable Hamiltonain. The latter, a variant 
of Calogero-Sutherland (CS) Hamiltonian in two dimension, is a generator of the dynamics of N particles interacting 
via long-ranged two body and three body interactions and confined by a harmonic oscillator potential [Qj8) . A similar 
technique has already been applied to Hermitian operators which maps the eigenvalue distribution to a state of the 
CS Hamiltonian in one dimension the known particle correlators for the latter case are then used to determine 
the eigenvalue correlations. A detailed knowledge of the non-stationary states of the 2-dimensional Calogero system 
can therefore be useful in dealing with a variety of NHRM. 

As well-known, the CS Hamiltonian is a fully integrable system with particles evolving in an ordered way with 
respect to time JlO| ; this implies a strong correlation between various particle states at different times. Our mapping 
thus reveals a very interesting feature of the eigenvalues of the operators associated with both conservative and 
dissipative dynamics. The eigenvalues evolve in a highly ordered, correlated way as the degree or the nature of 
the complexity changes. This implies a connection between the statistical nature of the eigenvalues of two different 
complex systems. As the nature of the interaction in the corresponding CS system is 1/r 2 type, 1-dimcnsional for 
conservative and 2-dimensional for dissipative systems, regardless of the nature of the complexity, a great deal of 
universality among the physical properties (those related to eigenvalue correlations) of complex systems seems to 
be present. The universality in the eigenvalue statistics of the operators in the regime of weak non-Hcrmiticity was 
indicated by another study too . 

We consider an ensemble of N x N non-Hermitian matrices H defined by a Gaussian measure p{H ) where p(H, y) cx 

exp[— X)f=i St i(yki;sHki-s + x ki;sHki; S Hik; S )] = C 'p(H) with C as the normalization constant, y and x related to 
variances and covariances of various matrix elements. Here the subscript s on a variable refers to one of its components, 
i.e real (s = 1) or imaginary (s = 2) part, with (3 as total number of the components. The above choice of p is made 
so as to include a large class of NHRM ensembles (for example yu = N/(l — r 2 ), xm- s = %ik;s = (—1) s ~ 1 tN/ (I — r 2 ) 
give GUE (t = I), Ginibrc (r = 0) and the ensemble of complex anti-symmetric matrices, referred as GASE later on, 
r = — I |l2); see |ll| for the eigenvalue statistics for the cases with < r < 1). 

A non-Hermitian matrix can be diagonalized by a transformation of the type A = U HV with A as the matrix of 
eigenvalues Xj and U and V as the left and right eigenvector matrices respectively. Let us first consider the case of 
an ensemble of non-Hermitian complex matrices {(3 — 2). Here the eigenvalues Xj = X)r=i W~ 1 '\7>> m g enera l ; are 
distributed over an area in the complex plan. Let P(z,y) be the probability of finding eigenvalues Xi of H between 
Zi and Zi + dzi at a given y, 
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with P as the unnormalized distribution: P(z, y) = J f(z, z*)p(H, y)dH with f(z, z*) = Ilili $( z i ~ \)$( Z Z ~ K)- 

The degree of difficulty associated with solving the integral eq.(l) motivates us to seek another route. We attempt to 
obtain an evolution equation for P(z, y) for the eigenvalues moving on the complex plane due to changing distribution 
parameters. As in the Hermitian case ||, we consider a combination of the rates of change of P in the parametric 

space, namely the sum S, S = Ef=i Eu (7 + (-l) Sx kks) &ta;s af^r + Xkl-,a sfrrr > and attempt to express it in 



ykl ' s dy kl;s + Xkl > s dx M;h a 

terms of the rates of changes of P in the eigenvalue space. (The reason behind the choice of sum S is that it can be 
reduced to the Schrodinger equation for CS Hamiltonian). This would require a knowledge of the rates of change of 
the eigenvalues as well as the eigenvectors due to a small change in the matrix element H k i which can be given as 
follows, 
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The parametric-dependence of P in sum S enters only through p(H) and as 
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the sum S can be written 
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where hi- 



EUElxalb I fm^ H M;s P dH, d = |Ef=iE M (7 + (-l) S ^) an d r referring to the 



components of the eigencvalues. With the help of eq.(5), the first term on the right hand side of eq.(8) can further be 
simplified, Es=i Efc i hl\s = E n r ~BI~ ( z nrP)- By using eqs.(6,7), the second term can also be rewritten as follows: 
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where G = E,(-1) S E M 
eq.(8) will now give the sum Si 
eigenvalue derivatives of P, 
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and Atv(z) = Y[f <k (zj ~ z k)- A substitution of eq.(9) in 
S + G + C\P, a combination of various parametric derivatives, in terms of the 
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However the sum S\ , on the other hand, can also be expressed as a derivative of P with respect to a single parameter 
Y = Y(y k i- S , Xki-s) where Y is a function of all y k i- s and x k i- s , given by the condition that 



Si 





EE 

s=l fc,Z 



Akl;. 



dP 

dyik;s 



Bhiis 



dP 



' dx 



Id:. 



C X P 



dP 
&Y 



•CiP 



1 dPi 
Ch~dY 



(11) 



2 



where A kl . s = yki-Al+Zi-lYxik^} and B u . s = [-yx k i. s + (-l) s x k i. s xi k . s + (~l) s y k i. s yi k . s ], C 2 = el CldY and Pi = C 2 P. 
The form of Y, fulfilling the desired condition, can therefore be obtained by solving following equations || (for all k, I 
and s values): 
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which gives Y = (1/N 2 ) , 5Zf =1 F(ykl;s) + with Y given by the initial conditions. Here F(y k i- S ) 
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-l) s Cfci ;s ) and constants 
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c fci;s and c fei;s given by relations y ik . s = c k i- s y k i- s and x 2 kl . s + (-l) s jx kl . s - c fc ; ;s y^. s - {-l) s c k i- s ykl;s = 0. 

. Here various y k i- s being indicators of the complexity of the system, Y can be termed as a complexity parameter. 
Now, by comparing the two forms of S% , the evolution of eigenvalues in terms of the parameter Y can be obtained 
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with (3 — 2 and Pi is related to the normalized distribution by P = CP1/C2. Note the analogy of the above 
equation to that of Hermitian case B but the evolution is now occuring on a complex plane. 



The steady state of eq.(13), P s = \Q 



n\ 2 — Ylj <k |Ajv(z)| 2 e 2 Sfe' 2fc ' , corresponds to |^ — > which can occur 
when almost all yki- s — > ^V/(l — r 2 ) and almost all Xfc; ;s — > (— 1) s ./Vt/(1 — r 2 ) with r — > 0, ±1. Here each r value 
leads to a different steady state, namely, GBE (r = 0) GUE (r = 1) and GASE (r = —1) with the distribution P s 
representing all the three cases. Note P s in each case agrees well with the nature of the matrix H for these limits 
which is complex for r = 0, complex hermitian for t = 1 (therefore real eigenvalues) and complex antisymmetric for 
r = — 1 (thus eigenvalues in equal and opposite pairs). 

The eq.(13) describeis a transition from a given initial ensemble (with Y — Y ) to either GBE, GUE or GASE 
with Y — Yq as the transition parameter. The nonequilibrium states of these transitions, given by non-zero finite 
values of Y — Yq, are various ensembles of the complex matrices corresponding to varying values of t/m's and a:;y's 
thus modelling different complex systems. Note the eq.(13) for Pi = Pi(/j,,Y\iio,Yq) has been obtained for arbitrary 
initial conditions, say Pi(fJ>o, Y ); the eigenvalue distribution Pi(/U, Y) = J Pi(/x, Y\fi , Y))Pi{h L o, Yo) of a given RNHE 
can therefore be found by solving the eq.(13) by using a convenient initial ensemble. Just as in the Hermitian case, 
the " convenience" depends on mathematical tractability of the integrals as well as on involved physics || . 

The case of Non- Hermitian real matrices ((3 = 1) can similarly be treated. Here eigenvalues are either real 
or form complex conjugate pairs and therefore if U n is an eigenvector corresponding to the complex eigenvalue 
A„, its complex conjugate will correspond to an eigenvector U*. Consider the case with L real and M complex 
conjugate pairs of the eigenvalues with N = L + 2M. The rates of change of the eigenvalues and the eigen- 
vectors are still given by eqs. (2-7) with H k ia = H k i. The distribution P in this case can be described by 

p = j* Iljli f( z ,z*)g(z,z*)p(H)dH with f{z,z*) = n^i^Mj - Zj)S { fj,j - z*), g(z,z*) = U^l+i 6 ^ ~ zj) 5 ^*] ~ 
Zj + M)5(^j+M — z j)Hp*j+M ~ Z j+Ai)- (A s obvious, here first L eigenvalues are chosen to be real and rest of them com- 
plex conjugate). Proceeding similarly as for the complex case, using eqs. (2-7) and equalities = — Y^nt^_ M ^ z "^ 
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where Y is still given by the eq.(ll) with (3 = 1 {y k i-i = yu and x a -i = x k i); Y = J^k.i F (Vki) with F(y k i) 
J dy k i{y k iVW)- 1 =)n 



and x 2 kl - jxu 
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, W = 7 2 + 4j/fej(cfcjJ/fcj - hi) and c ki and c fe z given by relations y lk = c kl y k i 



The steady state again occurs for r 
1/2 



0, ±1 with solution of eq.(14) is 



z*) 



The distribution P(z;t = 0,±1) is in agreement with the results 
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obtained in 113] by a different method. 

The eq.(13) for P(fi,Y) can be used to obtain n th order density correlator R n (z\, ..z n ; Y), defined by R n = 
(jv-n)i / P{z 1 Y)dz n+ i.AzM with dz n = dz„idz„2- The similar forms of the equations for P, in Hermitian and 
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non-Hermitian case, result in the same for the equations for R n too |9| . Again a direct integration of F-P equation 
(13) leads to the hierarchic relations among R n 
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For real applications, it is important to consider the limit N — > oo for fixed n. For p(z) ~ N 1 Ri(z) which fixes 
the scale for the eigenvalue fluctuations, the large N limit of eq.(15), gives the following form (with z = Ne) 
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where P refers to the priniciple part of the integral. The eq.(16) is valid under the similar approximations as in the 
Hermitian case, that is, by neglecting contributions from the terms containing 2 nd order cluster functions (see page 
142 of also |l^|) and the diffusion term; both are of the order of N 2 or lower. 

For n > 1, the correlators should be unfolded (that is, a rescaling of eigenvalues to result in a unit mean spacing) 

as follows: i?„(Ci, ..,Cni A) = LimiV — > oo r J^y)'' , rI{z' > -Y) wrtn C = / R\(z:Y)dz. Then, for large N and fixed Y, 
the left hand side of eq.(15) is 0(i?" ) while the right hand side is 0(i?™ +2 ). The transition therefore takes place for 
finite values of YR\ and a smooth transition can only be seen in terms of a parameter A = (Y — Yq)/D 2 (D = R± , 
the mean level spacing) pq| . Assuming that for Y — Yq, p is not singular (nor zero) and the i? n are well defined, 
the p, given by eq.(16), remains nearly unchanged for finite A. Keeping only 0(i?™ +2 ) terms, the eq.(15) can then be 
reduced to following form (see page 145 of [pill) 
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here b — 2 (for simplification, 7 is chosen to be unity). The above equation is obtained from eq.(15) by neglecting 
the linear drift of the eigenvalues which is dominated, by a factor N, by their diffusion and mutual repulsion. In fact 
the linear restoring force, responsible for the global behaviour p of the density of levels is entirely negligible on scales 
at which local fluctuations occur. On the other hand, the diffusion is ineffective on the global scale (see eq.(16)). As 
discussed above, the transition for R n occurs on the scales determined by Y « D 2 , while, for R\, the corresponding 
scale is given by Y « N D 2 . This indicates a clear separation of the scales of the global and local behaviour of the 
density. (It is worth noting here the similarity of eq.(17) to its Hermitian counterpart; however it does not imply the 
similarity of correlations, the former being on the complex plane). The hierarchical equation of correlations for the 
real assymetric case can be obtained by integrating eq.(14) which will again lead to a relation similar to eq.(17) but 
now (3=1, Yl r *° ^ e dropped and Q r replaced by Q. 

For n = 2 and small values of £, the integral term in eq.(17) makes a negligible contribution thus leading to following 
approximated closed form equation for R2(r = £2 — Ci) (with r = r\ + ir%) 
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which gives Rz(r) 



for small-r; for large r behaviour, it may be easier to consider the fourier transform of 



eq.(17). The hierarchic equation can then be used to obtain an approximate form of the higher order correlations. For 
example, the approximate information about R3 may be extracted by a substitution of large and small r behaviour 
of i?2 in eq.(15) with n = 2. 

An alternative route to obtain correlations is by exploiting the connection of eq.(13) to CS Hamiltonian. This 



can be shown by using the transformation $ = P/\Qn\^^ 2 in eq.(13) reducing it in a form §^ 
'Hamiltonian' H turns out to be the CS Hamiltonian in two dimensions (for simplification take 7 
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and rki = \rki\- Here G = g (with g = 1 for NHRM case with all real eigenvalues and 
g = 2 for the complex NHRM) and, unlike the complex Hermitian case (G = 0, g = 1), the inverse square term does 
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not drop out for the complex non-Hermitian case. Further as Y — > oo, the particles are in their ground state ip with a 
distribution ip^ . As the mapping between the eigenvalues and paticles requires V'o = Ps which is possible if the particles 
are considered as bosons. The bosonic radial eigenstate and the eigenvalues of the Hamiltonian are well known, 
i'n = Il^cfc |ri-rj| A e~i2j r <=l 2 £ n w ithi„ as Laguerre's polynomial and energy E n = [An+N(N-1)A+2N]/2N,A = 

§. The "state" ip or P(p,Y\m h Y ) can then formally be expressed as a sum over the eigenvalues and 
eigenfunctions which on integration over the initial state P(ho,Yq) would lead to the joint probability distribution 
P(n,Y) and thereby static (at a single parameter value) density correlations R n . The above correspondence can 
also be used to map the multi-parametric correlations of levels to multi-time correlations of the particle-positions 
. Although the explicit calculations of correlations, involves technical handling of various integrals and is still an 
unfinished task, nontheless our study reveals an important connection. The level correlations of different complex 
systems need not be studied separately, a thorough probing of the particle correlations of CS will give all the required 
information. The CS system being integrable in nature, the semiclassical techniques can also be very successful for 
the probing. 

The reasons for the correspondence between Gaussian NHRM and CS Hamiltonian are worth paying attention. 
Note the analogy with harmonic oscillator type confining potential in the CS system results from the Gaussian nature 
of the ensemble. The correspondence with 1/r 2 term comes from the mutual repulsion between eigenvalues. Its 
mathematical origin lies in the transformation from matrix space to eigenvalue-eigenvector space which is same for 
all the non-hermitian ensembles (belonging to same symmetry class irrespective of matrix element distribution). It 
should be possible, therefore, to map the non-Gaussian NHRM also to a variant of CS Hamiltonian, although with a 
different type of confining potential. For p{H) cx e~^ H ^ , the correspondence can be shown by following the similar 
steps as used for the Hermitian case [||. 

In this paper, we have studied the statistical properties of the eigenvalues of non-Hermitian systems. We find 
that the evolution of the eigenvalues is governed by an equation in which system-dependence enters only through the 
evolution parameter Y — Yq related to complexity of the system. It is possible that widely different systems with 
different values of the distribution parameters may share a same Y — lo~ va hie. Such systems will thus have similar 
statistical featurs which indiactes an underlying universality in the distribution of the eigenvalues of non-hermitian 
operators. Furthermore as the eigenvalue distribution for each complex system appears as a general state of CS 
system, any two such states, for example tpO^i) an d ipO^), being related by "time" evolution operator U(Y2,Yi), the 
eigenvalue distributions of the complex systems corresponding to Y\ and Yq. will also be connected. This would also 
reflect in their physical properties based on spectral fluctuations e.g conductance (assuming existence of ergodicity, 
that is, ensemble averages same as spectral averages). 

A detailed investigation of CS hamiltonian in arbitrary dimension can therefore give a lot of useful information 
about variety of complex systems and is therefore very much desirable 
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